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Abstract 

Quantum inequality restrictions on the stress-energy tensor for negative energy 
are developed for three and four-dimensional static spacetimes. We derive a general 
inequality in terms of a sum of mode functions which constrains the magnitude and 
duration of negative energy seen by an observer at rest in a static spacetime. This 
inequality is evaluated explicitly for a minimally coupled scalar field in three and four- 
dimensional static Robertson- Walker universes. In the limit of vanishing curvature, 
the flat spacetime inequalities are recovered. More generally, these inequalities contain 
the effects of spacetime curvature. In the limit of short sampling times, they take the 
flat space form plus subdominant curvature-dependent corrections. 
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1 Introduction 



Recently, several approaches have been proposed to study the extent to which quantum 
fields may violate the "weak energy condition," (WEC) T^^u^u'^ > 0, for all timelike vectors 
m'^. Such violations occur in the Casimir effect and in quantum coherence effects, where the 
energy density in a region may be negative. Violations of the local energy conditions first led 
to the averaging of the energy condition over timelike (or null) geodesies [|l|. This eventually 
led to the "averaged weak energy condition" , 



Here the total energy density averaged along an entire geodesic is constrained to be positive. 
While this does preserve in some form the weak energy condition, it does not constrain the 
magnitude of the negative energy violations from becoming arbitrarily large over an arbitrary 
interval, so long as there is compensating positive energy elsewhere. In cases such as the 
Casimir effect, where there is a constant negative energy density present in the vacuum state, 
Eq. fails. However, it is often still possible to prove "difference inequalities" in which 
Tf^i, is replaced by the difference in expectation values between an arbitrary state and the 
vacuum state ^. 

Quantum inequality (QI) type relations have been proven ^, ^ which do constrain the 
magnitude and extent of negative energy. For example, in 4D Minkowski space the quantum 
inequality for free, quantized, massless scalar fields can be written in its covariant form as 



for all values of tq and where is a timelike vector. The expectation value () is with 
respect to an arbitrary state and Tq is the characteristic width of the sampling function, 
To/h{t'^ + Tq), whose time integral is unity. Such inequalities limit the magnitude of the 
negative energy violations and the time for which they are allowed to exist. In addition, 
such QI relations reduce to the usual AWEC type conditions in the infinite sampling time 
limit. 

Flat space quantum inequality-type relations of this form have since been applied to 
curved spacetimes |^ by keeping the sampling time shorter than a "characteristic" cur- 
vature radius of the geometry. Under such circumstances, it was argued that the spacetime 
is approximately fiat, and the inequalities could be applied over a small region. In the 
wormhole geometry this led to wormholes which were either on the order of a Planck 
length in size or with a great disparity in the length scales that characterize the wormhole. 
In the "warp drive" geometry ^ it was found that if one wanted superluminal bubbles of 
macroscopically useful size then the bubble wall thickness would be of the order of a Planck 
length 0]. Although the method of small sampling times is useful, it does not address the 
question of how the curvature would enter into the quantum inequalities for arbitrarily long 
sampling times. 

In this paper we will address this issue. We will begin with a generalized theory in 
Section 2 which will allow us to find quantum difference inequalities in a globally static, 
but arbitrarily curved spacetime Then in Section 3 we will look at the case of the S"^ x R 
spacetime, the 3D equivalent to the Einstein universe. We will show that the QI's are 




(1) 




(2) 
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modified by a "scale function" wliicli is dependent upon tlie ratio of sampling time to the 
curvature radius {to /a). In Section 4, we will proceed to find similar QI's for the three cases 
of the 4D static Robertson- Walker spacetimes. In conclusion, we will discuss how AWEC 
type conditions in these particular spacetimes can be found. We will use units in which 

n = c = G = i. 



2 General Theory 

In this section, we will develop a formalism to find the quantum inequalities for a massive, 
minimally coupled scalar field in a generalized spacetime. This method will be applicable in 
globally static spacetimes, allowing us to use a separation of variables of the wave equation, 
and write the positive frequency mode functions as 

/A(x,t) = f/,(x)e-'-*. (3) 

The label A represents the set of quantum numbers necessary to specify the state. Addition- 
ally, the mode functions should have unit Klein-Gordon norm 

{fx,fx')=Sxy. (4) 

The above separation of variables can be accomplished when dt is a timelike Killing 
vector. Such a spacetime can be described by a metric of the form 

ds"^ = —dt^ + gij{x.)dxidxj. (5) 

Here Qij is the metric of the spacelike hypersurfaces that are orthogonal to the Killing vector 
in the time direction. With this metric, the wave equation is 

□0 _ = _52^ + ^ (d,J\^\g'^d,4>) - /iV = 0, (6) 



where g = det{g^u). The scalar field can then be expanded in terms of creation and 
annihilation operators as 

= EK/a + 4/a*), (7) 

A 

when quantization is carried out over a finite box or universe. If the spacetime is itself 
infinite, then we replace the summation by an integral over all of the possible modes. The 
creation and annihilation operators satisfy the usual commutation relations 0. 

In principle, quantum inequalities can be found for any geodesic observer 0. In the 
Robertson- Walker universes, static observers see the universe as having maximal symmetry. 
Moving observers lose this symmetry, and in general the mode functions of the wave equation 
can become quite complicated. Thus, for this paper we will concern ourselves only with static 
observers, whose four-velocity is the timelike Killing vector. The energy density that such 
an observer samples is given by 
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Upon substitution of the above mode function expansion into Eq. (H), one finds 



Too = Re^ W[f/^f/ve+^( 



a^ay 



U.Uye 



-i(w+uj')t 



AA' 



+ [dWldjUye+'^''-'''^'a\ay + dWxdjUye-'^'^+'^'^'ay.ay 



V [ulu. 



A'C 



a\ay + UxUye 



axay 



} 



(9) 



The last term is the expectation value in the vacuum state, defined by a\\Q) = for 
all A, and is formally divergent. The vacuum energy density may be defined by a suitable 
regularization and renormalization procedure, as will be discussed in more detail in Section 3 
and 4. In general, however, it is not uniquely defined. This ambiguity may be side-stepped 
by concentrating attention upon the difference between the energy density in an arbitrary 
state and that in the vacuum state, as was done in Ref. 0. We will therefore concern 
ourselves primarily with 

:Too:=Too-(0|Too|0), (10) 

where |0) represents the Fock vacuum state defined by the global timelike Killing vector. 
We will average the energy density as in Eq. (0), and find that the averaged energy density 
difference is given by 



Ap ^ 



n 



{: Too ■)dt 
t' + tl 



AA' 



ye 



-|w-a;'|to/„t 



{alay) + dmxdjUye 



-(a;+w')to 



{axay) 



V [uiUye-\---'\"^{a\ay) + UxUye-^^+-'^''^ {a^ay)] } 
We are seeking a lower bound on this quantity. It has been shown [|T^ that 



ReE/(A)*/(A')e 



{alay)>ReJ2fiXrfiX 



Mp-('^W)to/^t 



{alay). 



(11) 



(12) 



AA' 



AA' 



Upon substitution of this into Eq. (^TJ) we have 



Ap > ReY,{ujuj' [UlUy {alay) -UxUy {axay) 

AA' 

+ [d'U;d,Uy{a\ay) + dmxd,Uy{axay) 

+fi' [UlUy{a{ay) + UxUy{axay)] } e-(-+'^')*° . 



(13) 



We may now apply the inequalities proven in the Appendix. For the first and third term 
of Eq. (ID, apply Eq. ( |A7|) with hx = uUxe'^^^^ and hx = fiUxe~'^*°, respectively. For the 
second term of Eq. (|13D , apply Eq. (|AlD with Aij = gij and h\ = d'Uxe~'^^°- The result is 



Ap>-W {^lUlUx + d^U; d,Ux + p'UlUx) e 



(14) 
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This inequality may be re-written using the equation satisfied by the spatial mode functions: 

W,Ux + {uj^-fi^)Ux = 0, (15) 

to obtain 

A/3 > - E (^A + ^V^V.) |f/Ape-2-*o . (16) 

Here V* is the covariant derivative operator in the t = constant hypersurfaces. Therefore, 
given any metric which admits a global timelike Killing vector, we can calculate the limita- 
tions on the negative energy densities once the solutions to the wave equation in that curved 
background are known. 

Note that although the local energy density may be more negative in a given quantum 
state than in the vacuum, the total energy difference integrated over all space is non-negative. 
This follows from the fact that the normal-ordered Hamiltonian, 

: H := : Too : d!^x = ^uJx a\ax , (17) 

•' A 

is a positive-definite operator, so that {: H :) > 0. 

In the following sections, we will apply the general energy density inequality, Eq. ([T6|) , 
to three and four- dimensional Robertson- Walker universes. 



3 Quantum Inequality in a 3D Closed Universe 

Let us consider the three-dimensional spacetime with a length element given by 

ds"^ = -dt^ + (de^ + sin^ Odcp'^ 



Here constant time slices of this universe are two spheres of radius a. The wave equation on 
this background with a coupling of strength ^ to the Ricci scalar R is 

□0 - (/i^ + ^i?)0 = 0. (19) 

For the metric Eq. ([T8|), this equation becomes 

1 „ , . 1 ./2 



a^smt^ a^sm b' ^ J 

which has the solutions 

fUt.e,^) = ^l=YUO.^)e-'^'' / = 0,1,2,- •• (21) 

-I <m<l, (22) 

with eigenfrequencies 

cj; = a-y/(/ + 1) + 2^ + /iV . (23) 

The Yim are the usual spherical harmonics with unit normalization. The coupling parameter 
^ here can be seen to contribute to the wave functions as a term of the same form as the 
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mass. In this paper, however, we will look only at minimal coupling = 0). The lower 
bound on the energy density is then given by 



oo +1 



oo +1 



1=0 m=-l 



1=0 m=-l 



(24) 



However the spherical harmonics obey a sum rule |11 



+1 



21 



m=—l 



4:71 



(25) 



We immediately see that the difference inequality is independent of position, as expected 
from the spatial homogeneity, and that the second term of Eq. (plj) does not contribute. We 
have 

1 oo 

Ap> -^^(2/ + l)^;e"2-'*«. (26) 



1=0 



This summation is finite due to the exponentially decaying term. We are now left with 
evaluating the sum for a particular set of values for the mass /i, the radius a, and the 
sampling time to- Let 1] = tQ/a. Then 



Ap > 



2r/=^^(2/ + l)cu,e-2''^' 
1=1 



(27) 



where 



n{l + 1) + aV (28) 

The coefficient —1/(1671^11) is the right hand side of the inequality for the case of a massless 
field in an infinite 3-D Minkowski space. The "scale function", F{ri,fi), represents how the 
mass and the curvature of the closed spacetime affects the difference inequality. 



3.1 Massless Case 

In terms of the variable t] = to/ a, which is the ratio of the sampling time to the radius of 
the universe, we can write the above expression for F{ri, 0) as 

oo 

F{ri) = F{ri, 0) = 2ri'' ^(2/ + l)y7(7+l) e-W'('+i) (29) 

1=1 

A plot of F{ri) is shown in Figure 1. In the limit of r/ ^ 0, when the sampling time is very 
small or the radius of the universe has become so large that it approximates flat space, then 
the function F{ri) approaches 1, yielding the flat space inequality. . 

We can look at the inequality in the two asymptotic regimes of rj. In the large rj regime 
each term of greater I in the exponent will decay away faster than the previous. The / = 1 
term yields a good approximation. In the other regime, when the sampling time is small 
compared to the radius, we can use the Plana summation formula to calculate the summation 
explicitly: 

£ fin) + IfiO) = f f{x)dx + I f ^^^^^^^^^dx , (30) 
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where 



fix) = {2x + l)y^x(x + l)e-2V^(^+i). (31) 

Immediately we see that for our summation /(O) = 0. The first integral can be done with 
relative ease yielding 

r fix)dx = r(2x + + l)e-W^(^+i)rfx = -X, . (32) 

Jo Jo ^ 2r]-^ 

This term reproduces the fiat space inequality. The second integral in Eq. (^) therefore 
contains all the corrections due to non-zero curvature of the spacetime. Since 1] is small, we 
can Taylor expand the exponent around rj = 0, and then keep the lowest order terms. One 
finds 

f{ix) — f{—ix) 

dx ~ -Jo - r^/i + ■ ■ ■ , (33) 



where 



and 



(2x\/Vx^ + 1- X + \/Vx^ + l 



X 

V2^ ^ ^ dx ^ 0.265096 (34) 



'^-L (35) 



4x(2a;2 - 1) , 2 
— dx = 

Therefore the function F{ri) in the small r] limit is given by 

4 

F{r]) ~ 1 - 0.530192 r]^ + — r/^ + 0(r/^) + • • • (36) 

lo 

and in the large t] limit by 

F(r7) ^6V2r/3g-2v^r,^... _ ^g^^ 

Both of these asymptotic forms are plotted along with the exact form of F{rj) in Figure 
1. The graph shows that the two asymptotic limits follow the exact graph to a very large 
precision except in the interval of 1 < ?7 < 2. These results for the function F{ri), combined 
with Eq. (13), yield the difference inequality for the massless scalar field. 

The difference inequality does not require knowledge of the actual value of the renormal- 
ized vacuum energy, po- However, if we wish to obtain a bound on the energy density itself, 
we must combine the difference inequality and po- One procedure for computing po is anal- 
ogous to that used to find the Casimir energy in fiat spacetime with boundaries: one defines 
a regularized energy density, subtracts the the corresponding fiat space energy density, and 
then takes the limit in which the regulator is removed. A possible choice of regulator is to 
insert a cutoff function, g{uj), in the mode sum and define the regularized energy density as 

1 I oo 

Preg = -^-o Yl ^(^a) = ^— ^ Y.(2l + 1) Ul g{uji) . (38) 

In the limit that a — > oo, we may replace the sum by an integral and obtain the regularized 
flat space energy density: 

1 



PFSreg = ^ duj UJ g{uj) . (39) 
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F(ri) 




1 2 3 4 5 



Figure 1: 

Plot of the Scale Function F{ri) and its asymptotic forms. The solid hue is the exact form 
of the function. The dotted line is the small r) approximation, while the dot-dash line is the 
large rj approximation. 
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The renormalized vacuum energy density may then be defined as 



PO = lini(p,.eg - pFSreg) ■ (40) 

The vacuum energy density so obtained will be denoted by pcasimir- An analogous procedure 
was used in Ref. to obtain the vacuum energy density for the conformal scalar field in 
the 4D Einstein universe. An explicit calculation for the present case, again using the Plana 
summation formula, yields 

/o _ 0.265096 



This agrees with the result obtained by Elizalde using the zeta function technique. In 
this case pcasimir < 0, whereas the analogous calculation for the conformal scalar field in this 
three-dimensional spacetime yields a positive vacuum energy density, pcasimir = l/967ra^. 
It should be noted that this procedure works in this case only because the divergent part 
of preg is independent of a. More generally, there may be curvature- dependent divergences 
which must also be removed. 

Let us now return to the explicit forms of the difference inequality. In the limit that 
to > a, Eqs. and (|7|) yield 



3^2 



Ap > -^^e-2^*«/". (42) 



Similarly, in the limit that to ^ cl^ we find 

1 

lovrtn oOvra' 



- ~ 1 g_^3 ~" Pcasimir " H • (43) 



P-n>-T^-7^ + ---- (44) 



Thus the bound on the renormalized energy density in an arbitrary quantum state in the 
latter limit becomes 

1 

167rt|] 607ra^ 

From either of Eqs. (^31) or (^11) in the a ^ oo limit we obtain the fiat space limit. Here 
Po = so Ap = p. In 3D fiat space we can write the quantum inequality in a more covariant 
form as 



P = - _9 , _2 > -77— 3' (45) 



To {T^,U^'u'')dT ^ 1_ 

TT J-oo t"^ + To ~ IGtITq ' 

In the other limit when the sampling time becomes long, we find that Ap decays exponentially 
as a function of the sampling time. This simply refiects the fact that the difference in energy 
density between an arbitrary state and the vacuum state satisfies the averaged weak energy 
condition, Eq. (0). 

4 4D Robertson- Walker Universes 

Now we will apply the same method to the case of the three homogeneous and isotropic 
universes given by the four-dimensional static Robertson- Walker metrics. Here we have 

[e = 0] : Qijdx'dx^ = dx^ + dy^ + dz^ (46) 
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for the flat universe with no curvature (Minkowski spacetime). For the closed universe with 
constant radius a, i.e. the universe of constant positive curvature, the spatial length element 
is given by. 



[e = 1] : Qijdx^dx^ 



dx + sin^ x{dG + sin^ 



(47) 



where 0<x<^) < 6 < n, and < < 27r. The open universe [e = —1] is given by 
making the replacement sin x sinh x in Eq. (|47|) and now allowing x to take on the values 
< X < oo- To find the lower bound of the quantum inequalities above we must solve for 
the eigenfunctions of the covariant Helmholtz equation 



(48) 



where /i is the mass of the scalar field and uj\ is the energy. A useful form of the solutions 



for this case is given by Parker and Fulling WA . 



4.1 Flat and Open Universes in 4D 

One finds that in the notation developed in Section 2, the spatial portion of the wave func- 
tions for flat (Euclidean) space is given by (continuum normalization) 

[e = 0] : f/k(x) = [2cu(27r)3]-i/2eik-x^ (49) 



a; = ^|k|2 + ^2 ^50) 
k = (/ci, k2, ks) (—00 < kj < 00). 

It is evident that will be independent of position. This immediately removes the second 
term of the inequality in Eq. ([T6|) . 

In the open universe the spatial functions are given by 

[e = -l]: Uxi^) = {2a'u,)-'/'nl^j\x)Yim{e,^), (51) 



A = {qj,m). 

Here, < q < 00; / = 0, 1, ■ ■ ■ ; and m = —I, — Z + 1, ■ ■ ■ , +/. The sum over all states involves 
an integral over the radial momentum q. The functions Il^^\x) are given in Equation 5.23 
of Birrell and DaviesQ. Apart from the normalization factor, they are 

K^\x) oc sinh' x K^^^^^ j cosgx- (53) 
As with the mode functions of the 3-dimensional closed spacetime above, the mode functions 



of the open 4-dimensional universe satisfy an addition theorem |1J, |I5|, |16 

,2 



J:\W(x)YiUO,^)\' = ^. (54) 



Since the addition theorem removes any spatial dependence, we again get no contribution 
from the Laplacian term of the quantum inequality, Eq. ([T6|). Upon substitution of the mode 
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functions for both the flat and open universes into the quantum inequahty and using the 
addition theorem in the open spacetime case we have 



1 /""^ 

1^ = °!^ ^"s-Ii^sL*'*'^^^"""" 



(55) 



and 



1 

[e = -l]: Ap>-^y^ dqq'u^.e 



-2u) qtfj 



4,2,3 ,„ , (56) 

respectively. The 3-dimensional integral in momentum space can be carried out by making 
a change to spherical momentum coordinates. The two cases can be written compactly as 



Ap > 



(57) 



where 



uj{k, a, /i) = (ji"^ — e/a^ + /i^^ 



-1/2 



(5^ 



Note that e = for flat space and k = q/a for the open universe. This integral can be carried 
out explicitly in terms of modified Bessel functions Kn{z). The result is 



Ap > 



327r2t4 



z''K^{z)-z^K2{z] 



G{z) 



where 



2to\/p 



,2- 



(59) 



(60) 



The coefficient — 3/(327r^tQ) is the lower bound on p found in Refs. 0, |T^ for a massless 



scalar field in Minkowski spacetime. The function G{z) is the "scale function", similar to 
that found in the case of the 3D closed universe. It is the same function found in Ref. ||10| 
for Minkowski spacetime, (e = 0) and is plotted in Figure 2. Again we see that in the limit 
of z — > the scale function approaches unity, returning the fiat space massless inequality in 
four dimensions P. 



4.2 The Closed Universe 

In the case of the closed universe the spatial mode functions are the 4-dimensional spherical 
harmonics, which have the form [|14|, ^ 

[e = 1] : t/,(x) = (2c.„a=')-V2nW(x) Y,„(^,¥,), (61) 

A = (ra, /, m). 

Here, n = 0, 1, 2, ■ ■ ■ ; / = 0, 1, ■ ■ ■ , n; and m = — /, — / + 1, ■ ■ ■ , +/. The function 11*^"'') is 
found from n*^~) by replacing x by —ix Q by —i{n + 1) |l^. Alternatively, they can be 
written in terms of Gegenbauer polynomials [jl^ |18| as 

nit)(x)ocsin'xCi+_^,(cosx). (63) 
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G(z) 




z 

Figure 2: 

Plot of the Scale Function G{z) for the Open and Flat Universes. 
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In either case, the addition theorem is found from Eq. 11.4(3) of Ref.|[T8[] (with p = 2 and 
^ = rj). This reduces to 

E|nit^(x)Y,„(^,^)P = ^^^, (64) 

Im 

from which it is easy to show that the energy density inequahty, Eq. (p!6|), becomes 

1 



Ap > 



^(n+l)2cu„e-2-"*°. 



(65) 



n=0 



If we use the variable t] = to/a in the above equation, we can simphfy it to 

3 



Ap > 



i/(r/,p). 



(66) 



Again we find the flat space solution in 4 dimensions, multiplied by the scale function H{ri, p), 
which is defined as 



H{7], fi) = ^7]"" J^{n + l)yn(n + 2) + aV e-'V«("+2)+«V2^ 



(67) 



n=l 



and is plotted in Figure 3. The scale function here has a small bump in it, occurring at 
roughly t] ^ 0.5 with a height of 1.03245. This may permit the magnitude of the negative 
energy to be slightly greater for a massless scalar field in the Einstein universe than is allowed 
in a flat universe for comparable sampling times. A similar result was shown to exist for 
massive fields in a 2 dimensional Minkowski spacetime |TD|. 



4.3 Massless asymptotic limits in the Einstein universe 

As with the 3-dimensional closed universe, we can find the asymptotic limits of this function. 
We again follow the method of the previous section, assuming the scalar field is massless, 
and making use of the Plana summation formula to find 



where 



and 



h = {x + 1) Vx(x + 2) e" W-(-+2) 



2Re 



{ix + l)yia;(ix + 2)e-W^^(^^'+2) 



dx 



\x 



l)\J^/xF+A + x- 2x\J + 4 



X 



dx 



-0.356109 



(68) 
(69) 



(70) 
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The first integral can be done in terms of Struve, H„(2;), and Neumann, Nn{z), functions, 
with the resuh 



16 dt]'^ 



4 

H V . 

15 ' 



(71) 



If we follow the same procedure as in the previous section for defining the renormalized 
vacuum energy density for the minimally coupled scalar field in the Einstein universe, then 
we obtain 

_ 1 0.356109 

The same method yields pcasimir = l/4807r^a^ for the massless conformal scalar field ||12|| . 
Here our result for the minimal field, Eq. ([f2|), differs from that obtained by Elizalde ||13[ 
using the zeta function method, p'casimir = — 0.411502/47r^a'^. This discrepancy probably 
reflects the fact that the renormalized vacuum energy density is not uniquely deflned. The 
renormalized stress tensor in a curved spacetime is only deflned up to additional flnite renor- 
malizations of the form of those required to remove the inflnities. In general this includes the 
geometrical tensors ^^^H^i^ and ^"^^H^i^. (See any of the references in |jl9[ for the deflnitions 
of these tensors and a discussion of their role in renormalization.) In the Einstein universe, 
both of these tensors are nonzero and are proportional to 1/a'^. Thus the addition of these 
tensors to (T^j/) will change the numerical coefficient in pcasimir- The logarithmically di- 
vergent parts of {T^y) which are proportional to ^^^H^,^ and ^"^"'H^y happen to vanish in the 
Einstein universe, but not in a more general spacetime. In principle, one should imagine 
that the renormalization procedure is performed in an arbitrary spacetime, and only later 
does one specialize to a specific metric. (Unfortunately, it is computationally impossible to 
do this explicitly.) Thus the fact that a particular divergent term happens to vanish in a 
particular spacetime does not preclude the presence of finite terms of the same form. 

In the small 7] limit we can Taylor expand the Struve and Neumann functions in Eq. ([7l|) , 
and obtain 



+ ^(-0-356109)77^ + ^r]^ log r/ + 0(r/^) + 



(73) 



where 7 is Euler's constant, which arises in the Taylor series expansion of the Neumann 
function. This is similar to that of the 3-dimensional universe and again contains a term 
of the form of the Casimir energy. In the large rj limit we again keep just the first term of 
the series. Both of the asymptotic forms are plotted with the exact solution in Figure 3. 
We see that the asymptotic form is again a very good approximation except in the interval 
1 < < 2, as was the case for the 3-dimensional closed universe. 
The difference inequalities for massless fields are then given by 



Ap > 



2 , (3 + 47) 4^14, ^ 

7] + -T] log ?7 + 
O 



12 



PCa 



(74) 



for to <ti a and 



Ap > 



73 _ 



2V3to/a 



for tn ^ a . 



(75) 



Using Ap = Pren— PCasimir, whcre pcasimir is the expectation value of pren in the vacuum state, 
one can calculate the "total," or more formally the renormalized energy density that would 
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H(il) 




~ 1 2 3 4 

Figure 3: 

The Scale Function H{r]) for the Closed Universe. The solid line is the exact result. The 
dotted line is the asymptotic expansion for small r], while the dot-dash curve is the large rj 
approximation. The maximum occurs at ~ 0.5. 
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be constrained by the quantum inequalities, subject to renormalization ambiguities. For 
example in the 4-dimensional Einstein universe the renormalized energy density inequality 
is constrained by 



Pr 



> 



1 9 

1+3.^+ 



(3 + 4 7), 

12 



(76) 



for to ^ o-- Here the coefficient of the rj* term would be altered by a finite renormalization, 
but the rest of the expression is unambiguous. Similar expressions could be found for any of 
the other cases above. In the case of the fiat Robert son- Walker universe, there is no Casimir 
vacuum energy. Under such circumstances the difference inequality and the renormalized 
energy density inequality are the same, and are free of ambiguities. 



5 Discussion 

In the previous sections, we have derived a general form for a difference inequality for a 
massive, minimally coupled scalar field in a static curved spacetime. This general inequality 
was then evaluated explicitly in a 3-dimensional closed universe and in 4-dimensional closed 
and open universes. The resulting inequalities include the effects of the spacetime curvature. 
In the small sampling time limit, they reduce to the fiat space forms plus curvature-dependent 
corrections. These results lend support to the arguments given in Ref. to the effect that 
this should always be the case in curved spacetime. 

An additional feature of these difference inequalities is that they lead to averaged weak 
energy type conditions in the infinite sampling time limit. We found that in general 

to r ^00 ■)dt ^ a. 



^ 71 J-oc t^+tl - ^ ^ ' 

where n is the dimensionality of the spacetime, a„ is the coefficient of the fiat space inequality 
in n dimensions, and Fn{fi) is the scale function. In all of the cases above the scale function 
decays exponentially to zero in the large t] limit, which is the same as the large to limit. In 
the to ^ 00 limit the above expression reduces to 

/OO /"OO 
(: Too ■.)dt = / mTo,\ij) - (0|Too|0)] dt > 0, (78) 
-00 J —00 

where is an arbitrary state. 

The difference inequalities are derived directly from quantum field theory and are in no 
way dependent upon the standard uncertainty relations. It therefore appears that quantum 
field theory itself leads to constraints on negative energy densities (fluxes) without any 
apriori assumptions. In addition, the quantum inequalities in both fiat space and these 
generalized curved spacetimes are more restrictive on the magnitude and duration of the 
allowable negative energy densities than is the averaged weak energy condition. Finally it 
appears that the quantum inequalities themselves lead directly to AWEC type conditions in 
many spacetimes. 
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Appendix 



In this appendix, we wish to prove the following inequality: Let Ay be a real, symmetric 
n X n matrix with non-negative eigenvalues. (For the purposes of this paper, we may take 
either n = 2, for 3D spacetime, or n = 3, for 4D spacetime.) Further let h\ be a complex 
ra- vector, which is also a function of the mode label A. Then in an arbitrary quantum state 
IV'), the inequality states that 



\,X' ^ A 

In order to prove this relation, we first note that 

n 



(Al) 



(A2) 



where the V^;'^^ are the eigenvectors of A^^ and the > are the corresponding eigenvalues. 
Now define the hermitian vector operator 



A 



(A3) 



Note that 



n n 



a=l 



a=l 



Furthermore, 



Q'^A.Q^) = 2ReJ2A, [h'x*hi{aiay) + h',hi,{axay)]+J2Ajh\*hi, (A5) 

A,A' A 

from which Eq. ( |A1| ) with the '-|-'-sign follows immediately. The form of Eq. ( [AT| ) with the 
'—'-sign can be obtained by letting hx — ^ ihy 

As a special case, we may take Aij = Sij and obtain 



Re^ [hi ■ hy{a{ax') ± hx ■ hyiaxay) > I^a^ • 



(A6) 



A,A' 



As a further special case, we may take the vector hA to have only one component, e.g. 
hA = {hx, 0, 0), in which case we obtain 



Re^ \jix*hy{alay) ± hxhy{axay) 



A,A' 



(A7) 



This last inequality was originally proven in Ref. ||20[ for real hx, and a simplified proof 
using the method adopted here is given in Ref. [p^ |. 
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